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We report the characterization of a two-qubit processor implemented with two capacitively coupled 
tunable superconducting qubits of the transmon type, each qubit having its own non-destructive 
single-shot readout. The fixed capacitive coupling yields the ^/iSWAP two-qubit gate for a suitable 
interaction time. We reconstruct by state tomography the coherent dynamics of the two-bit register 
as a function of the interaction time, observe a violation of the Bell inequality by 22 standard 
deviations after correcting readout errors, and measure by quantum process tomography a gate 
fidelity of 90%. 



Quantum information processing is one of the most 
appealing ideas for exploiting the resources of quantum 
physics and performing tasks beyond the reach of classi- 
cal machines |T]. Ideally, a quantum processor consists 
of an ensemble of highly coherent two-level systems, the 
qubits, that can be efficiently reset, that can follow any 
unitary evolution needed by an algorithm using a uni- 
versal set of single and two qubit gates, and that can be 
readout projectively. In the domain of electrical quantum 
circuits [5], important progress fS'-'T has been achieved 
recently with the operation of elementary quantum pro- 
cessors based on different superconducting qubits. Those 
based on transmon qubits |3l HI El |9] are well protected 
against decoherence but embed all the qubits in a single 
resonator used both for coupling them and for joint read- 
out. Consequently, individual readout of the qubits is not 
possible and the results of a calculation, as the Grover 
search algorithm demonstrated on two qubits [3J , cannot 
be obtained by running the algorithm only once. Further- 
more, the overhead for getting a result from such a pro- 
cessor without single-shot readout but with a larger num- 
ber of qubits overcomes the speed-up gain expected for 
any useful algorithm. The situation is different for pro- 
cessors based on phase qubits [SI El [TT] , where the qubits 
are more sensitive to decoherence but can be read indi- 
vidually with high fidelity, although destructively. This 
significant departure from the wished scheme can be cir- 
cumvented, when needed, since a destructive readout can 
be transformed into a non-destructive one at the cost of 
adding one ancilla qubit and one extra two-qubit gate for 
each qubit to be read projectively. Moreover, energy re- 
lease during a destructive readout can result in a sizeable 
cross-talk between the readout outcomes, which can also 
be solved at the expense of a more complex architecture 

unmi]. 

In this work, we operate a new architecture that comes 
closer to the ideal quantum processor design than the 
above-mentioned ones. Our circuit is based on frequency 
tunable transmons that are capacitively coupled. Al- 
though the coupling is fixed, the interaction is effective 
only when the qubits are on resonance, which yields the 
^/iSW AP universal gate for an adequate coupling dura- 



tion. Each qubit is equipped with its own non-destructive 
single-shot readout [T^ and the two qubits can be 
read with low cross-talk. In order to characterize the cir- 
cuit operation, we reconstruct the time evolution of the 
two-qubit register density matrix during the resonant and 
coherent exchange of a single quantum of excitation be- 
tween the qubits by quantum state tomography. Then, 
we prepare a Bell state with concurrence 0.85, measure 
the CHSH entanglement witness, and find a violation of 
the corresponding Bell inequality by 22 standard devia- 
tions. We then characterize the y/iSWAP universal gate 
operation by determining its process map with quantum 
process tomography [IJ. We find a gate fidelity of 90% 
due to qubit decoherence and systematic unitary errors. 

The circuit implemented is schematized in Fig|TjL: the 
coupled qubits with their respective control and read- 
out sub-circuits are fabricated on a Si chip (see sup- 
plementary information SI). The chip is cooled down 
to 20 mK in a dilution refrigerator and connected to 
room temperature sources and measurement devices by 
attenuated and filtered control lines and by two mea- 
surement lines equipped with cryogenic amplifiers. Each 
transmon j = /, // is a capacitively shunted SQUID 
characterized by its Coulomb energy for a Cooper 
pair, the asymmetry dj between its two Josephson junc- 
tions, and its total effective Josephson energy Ej{(l)j) — 
E^j |cos(a;j)| ^1 + d?ta.n^{xj), with xj — 7r0j/0o, 0o the 

flux quantum, and the magnetic flux through the 
SQUIDs induced by two local current lines with a 0.5 GHz 

bandwidth. The transition frequencies i/j ~ 2Ej-,E-j/h 
between the two lowest energy states |0)^ and can 
thus be tuned by . The qubits are coupled by a capac- 
itor with nominal value Cc — 0.13 fF and form a register 
with Hamiltonian H — h (— i^jcr^ — ^ii^^^} + '^9<^\<^^^ /2. 
Here h is the Planck constant, ctx, y,z are the Pauli opera- 
tors, 2g = \J E\,E^viVii/ Ecc ^ is the couphng fre- 
quency, and Eqc the Coulomb energy of a Cooper pair on 
the coupling capacitor. The two-qubit gate is deflned in 
the uncoupled basis {\uv)} = {|0)j , |l)j} ® {|0)jj , |l)ii}. 
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Figure 1: (a): circuit scliematics of the experiment with 
qubits in green and readout circuits in grayed blue, (b) Left 
panel: Spectroscopy of the sample showing the resonator fre- 
quencies v]^^ (horizontal lines), and the measured (disks, tri- 
angles) and fitted (lines) qubit frequencies viji as a function 
of their flux bias (piji when the other qubit is far detuned. 
Right panel: Spectroscopic anticrossing of the two qubits re- 
vealed by the 2D plot of poi + Pio as a function of the probe 
frequency and of ^i, at i^n — 5.124 GHz. (c) Typical pulse 
sequence including X or Y rotations, a V iSWAP" gate, Z 
rotations, and tomographic and readout pulses. Microwave 
pulses a{t) for qubit (green) and for readout (blue) are drawn 
on top of the i^i,ii(<^) dc pulses (red lines). 



at a working point Mi^u where the qubits are sufficiently 
detuned (fn — 3> 2g) to be negligibly coupled. Bring- 
ing them on resonance at a frequency v in a time much 
shorter than l/2g but much longer than l/i^, and keeping 
them on resonance during a time At, one implements an 

operation Qi.Qu.y/iSWAP ''^^ which is the product 
of the 
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gate to an adjustable power and of two single qubit phase 



gates 8j = exp (i6'j(jj/2) accounting for the dynamical 
phases 9j = J 27r(j/ — i'j)dt accumulated during the cou- 
pling. The exact \/iSWAP gate can thus be obtained 
by choosing At — l/8g and by applying a compensation 
rotation QJ'^ to each qubit afterward. 

For readout, each qubit is capacitively coupled to its 
own A/2 coplanar waveguide resonator with frequency 
ly^^ and quality factor Qj. This resonator is made non 
linear with a Josephson junction and is operated as a 
Josephson bifurcation amplifier, as explained in detail in 
[T?] . The homodyne measurement (see Fig|T^) of two mi- 
crowave pulses simultaneously applied to and refiected 
from the resonators yields a two-bit outcome uv that 
maps with a high fidelity the state \uv) on which the reg- 
ister is projected; the probabilities puv of the four possible 
outcomes are determined by repeating the same experi- 
mental sequence a few 10* times. Single qubit rotations 
u {9) by an angle 9 around an axis u of the XY plane 
of the Bloch sphere are obtained by applying Gaussian 
microwave pulses through the readout resonators, with 
frequencies i^j, phases ipj = (]l,lt), and calibrated area 
Aj oc 9. Rotations around Z are obtained by changing 
temporarily vi u with dc pulses on the current lines. 

The sample is first characterized by spectroscopy (see 
FigjljD) and a fit of the transmon model to the data yields 
the sample parameters (see S2). The working points 
where the qubits are manipulated (M^'^^), resonantly 
coupled (C), and read out {R^'^^) are chosen to yield suf- 
ficiently long relaxation times ^ 0.5 jis during gates, neg- 
ligible residual coupling during single qubit rotations and 
readout, and best possible fidelities at readout. Figure [TJa 
shows these points as well as the spectroscopic anticross- 
ing of the two qubits at point C, where 2g = 8.3 MHz 
in agreement with the design value of Cc- Then, read- 
out errors are characterized at R^'^^ (see Fig. S3.1): In 
a first approximation, the errors are independent for the 
two readouts and are of about 10% and 20% when read- 
ing |0) and |1) respectively. This limited fidelity results 
for a large part from energy relaxation of the qubits at 
readout. In addition we observe a small readout cross 
talk, i.e. a variation of up to 2% in the probability of 
an outcome of readout j depending on the state of the 
other qubit. All these effects are calibrated by measuring 
the four p^v probabilities for each of the four \uv) states, 
which allows us to calculate a 4 x 4 readout matrix TZ 
linking the Puv's to the \uv) populations. 

Repeating the pulse sequence shown in Fig|l]3 at — 
5.247 GHz, Af" = C = 5.125 GHz, R^ = 5.80 GHz, 
R^^ ~ 5.75 GHz, and applying the readout corrections 
TZ, we observe the coherent exchange of a single excita- 
tion initially stored in qubit /. We show in Fig|2] the 
time evolution of the measured \uv) populations, in fair 
agreement with a prediction obtained by integration of a 
simple time independent Liouville master equation of the 
system, involving the independently measured relaxation 
times TI — 436 ns and T^^ — 520 ns, and two indepen- 
dent effective pure dephasing times = = 2.0 /is 
as fitting parameters. Tomographic reconstruction of 
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Figure 2: Coherent swapping of a single excitation be- 
tween the qubits. (a) Experimental (solid lines) and fitted 
(dashed lines) occupation probabilities of the four computa- 
tional states |00) .. as a function of the coupling duration. 
No Z or tomographic pulses are applied here. (b,c) State to- 
mography of the initial state (left) and of the state produced 
by the V iSW AP gate (right), (b) Ideal (empty bars) and ex- 
perimental (color filling) expected values of the 15 Pauli op- 
erators XI, ..,ZZ. (c) Corresponding ideal (color filled black 
circles with pointer) and experimental (red circle and pointer) 
density matrices, as well as fidelity F and concurence C. Each 
matrix element is represented by a circle with an area propor- 
tional to its modulus (diameter = cell size for unit modulus) 
and a phase pointer giving its argument. 



the register density matrix p is obtained by measuring 
the expectation values of the 15 two-qubit PauH opera- 
tors {Pk} = {XI, .., ZZ}, the Xj and Yj measurements 
being obtained using tomographic pulses Ij^ (—90°) or 

jtj (90°) just before readout. The p matrix is calculated 
from the Pauli set by global minimization of the Hilbert- 
Schmidt distance between the possibly non-physical p 
and all physical (i.e. positive-semidefinite) p's. This can 
be done at regular interval of the coupling time to pro- 
duce a movie of p (At) (see supplementary on line ma- 
terial) showing the swapping of the |10) and |01) popu- 
lations at frequency 2g, the corresponding oscillation of 
the coherences, as well as the relaxation towards |00). 
Figure shows { (Pk) } and p only at At=0 ns and after 



iSWAP obtained at At=31 ns with ^rotations of 
and On ~ -1-60°. The fidelity F = {ipid\ p \i^id) 
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Figure 3: Test of the CHSH-Bell inequality on a |10)+e'''' jOl) 
state by measuring the qubits along X'or and XjJ or 
(see top-left inset), respectively. Blue (resp. red) error bars 
are the experimental CHSH entanglement witness determined 
from the raw (resp. readout errors corrected) measurements 
as a function of the angle ip between the measuring basis, 
whereas solid line is a fit using ip as the only fitting param- 
eter. Height of error bars is ± one standard deviation <j{N) 
(see bottom-right inset), with A'^ the number of sequences per 
point. Note that averaging beyond A'' = 10^ does not improve 
the violation because of a slow drift of ip. 



of p with the ideal density matrices \tpid) {tpid\ are 95% 
and 91%, respectively, and are limited by errors on the 
preparation pulse, statistical noise, and relaxation. 

To quantify in a different way our ability to entan- 
gle the two qubits, we prepare a Bell state |10) + e''^ |01) 
(with ^ — 6ii — 9i) using the pulse sequence of Fig|lJ; with 
At — 31 ns and no 0~^rotations, and measure the CHSH 
entanglement witness (XX^) + (XY^) + (YY^) - (YX^) 
as a function of the angle tp between the orthogonal mea- 
surement bases of qubit / and II. Figure [3] compares the 
results obtained with and without correcting the readout 
errors, with what is theoretically expected from the de- 
coherence parameters indicated previously: unlike in [lOJ 
and because of a readout contrast limited to 70 — 75%, 
the witness does not exceed the classical bound of 2 with- 
out correcting the readout errors. After correction, it 
reaches 2.43, in good agreement with the theoretical pre- 
diction (see also [2]), and exceeds the classical bound 
by up to 22 standard deviations when averaged over 10^ 
sequences. 

In a last experiment, we characterize the imperfec- 
tions of our V iSWAP gate by quantum process to- 
mography jl]. We build a completely positive map 
Pout = £{Pin) = J2m,nXmnPmPinPn^ characterized by a 
16 X 16 matrix x expressed here in the modified Pauli 

operator basis |Pk| = {I,X,Y' = iY,Z}'^'^, for which 

all matrices are real. For that purpose, we apply the 
gate (using pulse sequences similar to that of [l];, with 
At — 31 ns and 0~^rotations) to the sixteen input states 
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Figure 4: Map of the implemented \/ iSWAP gate yielding a 
fidelity of 90%. (a) Superposition of the ideal (empty thick 
bars) and experimental (color filled bars) lower part of the 
Hermitian matrix x (elements below 1% not shown). Each 
complex matrix element is represented by a bar with height 
proportional to its modulus and a red phase pointer at the 
top of the bar (as well as a filling color for experiment) giving 
its argument (top left inset). Expected peaks are marked by 
a star, (b) Lower part of the x error matrix (red circles - see 
text), with the same convention as in Fig[2j but with circles 
magnified for readability (a one-cell diameter represents a 8 
% modulus). Main visible contributions (continuous circles) 
are explained in text. 



{|0) , |1) , |0) + |1) , |0)+i and characterize both the 

input and output states by quantum state tomography. 
By operating as described previously, we would obtain 
apparent input and output density matrices including er- 
rors made in the state tomography itself, which we don't 
want to include in the gate map. Instead, we fit the 16 
experimental input Pauli sets by a model including am- 
plitude and phase errors for the X and Y preparation and 
tomographic pulses (see S4) , in order to determine which 
operator set {P^} is actually measured. The input and 
output matrices /9in,out corrected from the tomographic 
errors only are calculated by inverting the linear relation 
{{P§) = Tr (pPjJ)} and by applying it to the experimen- 
tal Pauli sets. We then calculate from the {pin.out} set an 
Hermitian x matrix that is not necessarily physical due 
to statistical errors, and which we render physical by tak- 
ing the nearest Hermitian positive matrix. This final x 
matrix is shown in Fig. |4]and compared to the ideal one, 
Xid, which yields a gate fidelity — Tr (x-Xid) — 0.9 
|15] . To better understand the imperfections, we also 
show the map x of the actual process preceded by the 
inverse ideal process |16] . The first diagonal element of x 
is equal to F by construction. Then, main visible errors 
arise from unitary operations and reduce fidelity by 1-2% 
(a fit yields a too long coupling time inducing a 95° swap 
instead of 90° and 81.2 rotations too small by 3.5° and 7° 
respectively). On the other hand the known relaxation 
and dephasing times reduces fidelity by 8% but is barely 
visible in x due to a spread over many matrix elements 
with modulus of the order of or below the 1 — 2% noise 
level. 

As a conclusion, we have demonstrated a high fidelity 
ViSWAP gate in a two josephson qubit circuit with in- 
dividual non-destructive single-shot readouts, observed a 
violation of the CHSH-Bell inequality, and followed the 
register's dynamics by tomography. Although quantum 
coherence and readout fidelity are still limited in this cir- 
cuit, they are sufficient to test in the near future simple 
quantum algorithms and get their result in a single run, 
which would demonstrate the concept of quantum speed- 
up. 
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technical support from P. Orfila, P. Senat, and J.C. Tack, 
as well as financial support from the European research 
contracts MIDAS and SOLID, from ANR Masquelspec 
and C'Nano and from the German Ministry of Education 
and Research. 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation [3] L. DiCarlo et al, Nature 460, 240 (2008). 

and Quantum Information (Cambridge University Press, [4] L. DiCarlo et al, Nature 467, 574 (2010). 

Cambridge, UK, 2000). [5] T. Yamamoto et al, Phys. Rev. B 82, 184515 (2010). 

[2] J. Clarke and F. Wilhelm, Nature 453, 1031 (2008). [6] R. C. Bialczak et al, Nature Physics 6, 409-413 (2010). 



5 



[7] J. M. Chow et al, Phys. Rev. Lett. 107, 080502 (2011). 

[8] J. Koch et al, Phys. Rev. A 76, 042319 (2007). 

[9] J. A. Schreier et al, Phys. Rev. B 77, 180502(R) (2008). 
[10] M. Ansmann et al, Nature 461, 504 (2009). 
[11] M. Mariantoni et al, Science, DOT: 10.1126/sci- 

ence.1208517 (2011). 
[12[ I. Siddiqi et al, Phys. Rev. Lett. 93, 207002 (2004). 
[13| F. Mallet et al, Nature Physics 5, 791 (2009). 
[14| J. M. Chow et al, Phys. Rev. A 81, 062325 (2010). 



[15] Note that Fg is also equal to Shumacher's fidelity 
Tr [Sj^s] /Tr [sISm] with S (resp. Su) the super oper- 
ator of the actual (resp. ideal unitary) process, and that 
fidelities F for the 16 outputs states range between 80% 
and 99.5% 

[16] A. G. Kofman and A. N. Korotkov, Phys. Rev. A, vol. 

81, no. 4, 042103 (2009), and private communication. 
[17[ F. R. Ong et al, Phys. Rev. Lett. 106, 167002 (2011). 



Supplementary material 

51. Sample preparation 

The sample is fabricated on a silicon chip oxidized over 50 nm. A 150 nm thick niobium layer is first deposited by 
magnetron sputtering and then dry-etched in a SFq plasma to pattern the readout resonators, the current lines for 
frequency tuning, and their ports. Finally, the transmon qubit, the coupling capacitance and the Josephson junctions 
of the resonators are fabricated by double-angle evaporation of aluminum through a shadow mask patterned by 
e-beam lithography. The first layer of aluminum is oxidized in a Ar — O2 mixture to form the oxide barrier of the 
junctions. The chip is glued with wax on a printed circuit board (PCB) and wire bonded to it. The PCB is then 
screwed in a copper box anchored to the cold plate of a dilution refrigerator. 

52. Sample parameters 

The sample is first characterized by spectroscopy (see Fig. l.b of main text). The incident power used is high 
enough to observe the resonator frequency i/r, the qubit line vqi, and the two-photon transition at frequency 1^02/2 
between the ground and second excited states of each transmon (data not shown). A fit of the transmon model 
to the data yields the sample parameters E]/h = 36.2 GHz, Ej^/h = 0.98 GHz, di ^ 0.2, Ef/h = 43.1GHz, 
E^(i/h = 0.87 GHz, du = 0.35, = 6.84 GHz, and = 6.70 GHz. The qubit-readout anticrossing at ly ^ 
yields the qubit-readout couplings ~ ~ 50 MHz. Independent measurements of the resonator dynamics (data 
not shown) yield quality factors Qi = Qu = 730 and Kerr non linearities [Hin] Ki/u^^ :^ Ku/vH ~ -2.3±0.5 x lO"'^. 

53. Experimental setup 

• Qubit microwave pulses: The qubit drive pulses are generated by two phase-locked microwave generators whose 
continuous wave outputs are fed to a pair of I/Q-mixers. The two IF inputs of each of these mixers are provided 
by a 4-Channell GS/s arbitrary waveform generator (AWG Tektronix AWG5014). Single-sideband mixing in 
the frequency range of 50-300 MHz is used to generate multi-tone drive pulses and to obtain a high ON/OFF 
ratio (> 50 dB) of the signal at the output of the mixers. Phase and amplitude errors of the mixers are corrected 
by measuring the signals at the output and applying sideband and carrier frequency dependent corrections in 
amplitude and offset to the IF input channels. 

• Flux Pulses: The flux control pulses are generated by a second AWG and sent to the chip through a transmission 
line, equipped with 40 dB of attenuation distributed over different temperature stages and a pair of 1 GHz 
absorptive low-pass filters at 4 K. The input signal of each fiux line is fed back to room temperature through 
an identical transmission line and measured to compensate the non-ideal frequency response of the line. 

• Readout Pulses: The pulses for the Josephson bifurcation amplifier (JBA) readouts are generated by mixing 
the continuous signals of a pair of microwave generators with IF pulses provided by a 1 GS/s arbitrary function 
generator. Each readout pulse consists of a measurement part with a rise time of 30 ns and a hold time of 100 
ns, followed by a 2/xs long latching part at 90 % of the pulse height. 

• Drive and Measurement Lines: The drive and readout microwave signals of each qubit are combined and sent to 
the sample through a pair of transmission lines that are attenuated by 70 dB over different temperature stages 
and filtered at 4K and 300 mK. A microwave circulator at 20 mK separates the input signals going to the chip 
from the reflected signals coming from the chip. The latter are amplifled by 36 dB at 4 K by two cryogenic 
HEMT ampliflers (GIT Gryo 1) with noise temperature 5K. The reflected readout pulses get further amplifled 
at room temperature and demodulated with the continuous signals of the readout microwave sources. The IQ 
quadratures of the demodulated signals are sampled at 1 GS/s by a 4-channel Data Acquisition system (Acqiris 
DG282). 
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S4. Readout characterization 



Errors in our readout scheme are discussed in detail in [13] for a single qubit. First, incorrect mapping |0) — > 1 
or|l) ^ of the projected state of the qubit to the dynamical state of the resonator can occur, due to the stochastic 
nature of the switching between the two dynamical states. As shown in FigjS] the probability p to obtain the outcome 
1 varies continuously from to 1 over a certain range of drive power Pd applied to the readout. When the shift 
in power between the two P|o>,|i> (^d) curves is not much larger than this range, the two curves overlap and errors 
are significant even at the optimal drive power where the difference in p is maximum. Second, even in the case of 
non overlapping P|o).|i) (-Pd) curves, the qubit initially projected in statejl) can relax down to |0) before the end of 
the measurement, yielding an outcome instead of 1. The probability of these two types of errors vary in opposite 
directions as a function of the frequency detuning A = — v > between the resonator and the qubit, so that a 
compromise has to be found for A. Besides, the contrast c = Max — P|o)) can be increased [12] by shelving state 
|1) into state |2) with a microwave tt pulse at frequency 1^12 just before the readout resonator pulse. The smallest 
errors Cq' and e{ when reading |0) and |1) are found for Aj = 440 MHz and An — 575 MHz and are shown by 
arrows in the top panels of Figlsl ef, — 5% and e\ — 13% (contrast c\ = \ — e\ — — 82%), and eg^ = 5.5% and 
e}!^ = 12% (cii = 82%). When using the |1) -> |2) shelving before readout, e[, = 2.5% and = 9.5% (contrast 
c\ == 1 — Bq — 63 = 88%), and ep — 3% and e}-^ = 8% (cn = 89%). These best results are very close to those obtained 
in [12[, but are unfortunately not relevant to this work. 

Indeed, when the two qubits are measured simultaneously, one has also to take into account a possible readout 
crosstalk, i.e. an influence of the projected state of each qubit on the outcome of the readout of the other qubit. We 
do observe such an effect and have to minimize it by increasing Ai.n up to ~ 1 GHz with respect to previous optimal 
values and by not using the shelving technique. An immediate consequence shown in Fig|5jb) is a reduction of the ciji 
contrasts. The errors when reading |0) and |1) are now ef, — 19% and e\ = 7% (contrast ci = 74%) and e^^ 19% 
and e^^ — 12% (contrast cn — 69%). Then to characterize the errors due to crosstalk, we measure the 4x4 readout 
matrix TZ linking the probabilities p^^ of the four possible uv outcomes to the population of the four \uv) states. As 
shown in Figjsj^c-d), we then rewrite TZ = Cqt- (Ci (E) Cu) as the product of a 4 x 4 pure crosstalk matrix Cqt with the 
tensorial product of the two 2x2 single qubit readout matrices 



/ 1 i,n i-ii \ 
'^1,11 — I,II , I. II 



We also illustrate on the figure the impact of the readout errors on our swapping experiment by comparing the 
bare readout outcomes uv, the outcomes corrected from the independent readout errors only, and the|uw) population 
calculated with the full correction including crosstalk. 

We now explain briefly the cause of the readout crosstalk in our processor. Unlike what was observed for other 
qubit readout schemes using switching detectors [5[ , the crosstalk we observe is not directly due to an electromagnetic 
perturbation induced by the switching of one detector that would help or prevent the switching of the other one. 
Indeed, when both qubits frequencies t'l n are moved far below the readout crosstalk disappears: the switching 
of a detector has no measurable effect on the switching of the other one. The crosstalk is actually due to the rather 
strong ac-Stark shift ~ 2 (nn — ni) g^/{R — v^i) ~ 500 MHz of the qubit frequency when a readout resonator switches 
from its low to high amplitude dynamical state with rtL ~ 10 and nn ~ 10^ photons, respectively. The small residual 
effective coupling between the qubits at readout can then slightly shift the frequency of the other resonatator, 
yielding a change of its switching probability by a few percent. Note that coupling the two qubits by a resonator 
rather than by a flxed capacitor would solve this problem. 



S5. Removing errors on tomographic pulses before calculating the gate process map 



Tomographic errors are removed from the process map of our ^JiSWAP gate using the following method. The 
measured Pauli sets corresponding to the sixteen input states are first fitted by a model including errors both in the 
preparation of the state (index prep) and in the tomographic pulses (index tomo). The errors included are angular 
errors efjj'^ on the nominal tt rotations around Xiji, 77}^jj'''*°™°and (5?™?^*°"^° on the nominal tt/2 rotations around 

Xiji and li,ii, a possible departure ^ui from orthogonality of (^Xi,}^i^ and (^Xii,Yii^, and a possible rotation fii^u 

of the tomographic XY frame with respect to the preparation one. The rotation operators used for preparing the 
states and doing their tomography are thus given by 
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Figure 5: Readout imperfections and their correction, (a) Switching probabihties of the readouts as a function of their driving 
power, with the qubit prepared in state|0) (blue), jl) ( red), or |2) (brown), at the optimal readout points. The arrows and 
dashed segments indicate the readout errors and contrast, at the power where the later is maximum, (b) Same as (a) but at 
readout points _R'''' used in this work, (c-d) Single readout matrices Ci.n and pure readout crosstalk matrix Cct characterizing 
the simultaneous readout of the two qubits. (e-g) bare readout outcomes uv, outcomes corrected from the independent readout 
errors only, andluz;) population calculated with the full correction including crosstalk for the swapping experiment of Fig|2] 
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Figure 6: Fitting of the pulse errors at state preparation and tomography. Measured (red) and fitted (blue - see text) Pauli 
sets (Pk) for the sixteen targeted input states {|0) , |f) , |0) + |1) , |0) + i |f)}®^. The {II , IX, lY, IZ, XI , ...} operators indicated 
in abscisse are the targeted operators and not those actually measured (due to tomographic errors). 
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L g-i(^/2+5rif )[cos(Ci,„)o-J,'"-sin(Ci.„)<Ti'"]/2^ 
j^tomo^^^2) = g-i(^/2+'Jij"°)[sin(/Ji,ii)cr^^'"+cos(/^iji)(Tj,-"]/2 
y^tomo^_^^2) — g+i('f/2+*l°i"°)[cos(/ii,ii+{i_n)(Tj,'"-sin(/ii,ii+4i,ii)c7j,'"]/2 

The sixteen input states are then {pf^ = U \0) {0\W} with {[/} = {/j, XjP"P(7r), yjP'''=P(7r/2), XjP'''P(-7r/2)} (g) 
{Jii,XP"P(7r),FP''°P(7r/2),XPj'''=P(-7r/2)}, and each input state yields a Pauh set {(P°) = Tr (p?„P^)} with {P^} = 
{h, X^, Y{, Zi} ® {/„, Xfi, Y[«i, Zii), X° = r'°-°(-7r/2)ta,r*°"-(-7r/2), and = X*™-(7r/2)ta,X*°"-(7r/2). Fig- 
urejG^shows the best fit of the modelled {{P^)} set to the measured input Pauli sets, yielding ffP'^'P ~ —1°, ePj'^^P = —3°, 
^prV^ 3=, j^P'-^P = 4°, ^P'-^P = -6°, 511°^ = -3°, ?7i'°™° = -6°, t];*?'"" = -4°, Af""^" = 12°, Ajp° = 5°, 6 = 1°, 6i = -2°, 
and /zi = /III = —11°. 

Knowing the tomographic errors and thus {(PjJ)}, we then invert the linear relation {{P^) — Tr (pP^)} to find the 

16 X 16 matrix B that links the vector {P^] to the columnized density matrix i.e. = B.{P^]. The matrix B is 
finally applied to the measured sixteen input and sixteen output Pauli sets to find the sixteen (pin,,Pout)k couples to 
be used for calculating the gate map. 



